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基于 8 面 截断 重要 抽样 法 可 靠 性 灵敏 度 估计 
及 其 方差 分 析 * 


Jk 峰 !2， 虽 震 宙 2， ÆRA? 

(1- 西北 工业 大 学 力学 与 土木 建筑 学 院 ， 西 安 710129; 2- 西北 工业 大 学 航空 学 院 ， 西 安 710072) 
摘 要 : 基于 6 面 的 截断 重要 抽样 法 可 以 用 来 求解 单 失效 模式 可 靠 性 灵敏 度 .该 方法 在 设计 点 处 作 失 效 重 
的 虚拟 切面 6 面 ， 而 8 面 将 变量 空间 分 割 成 重要 抽样 区 域 及 和 非 重 要 抽样 区 域 S， 在 及 和 SS 区 

域 分 别 建立 相应 的 截断 重要 抽样 密度 函数 hr(x) 和 hs(x)， 从 hr(x) 和 hs(x) 中 抽取 的 样本 量 

按照 及 和 SS 区 域 对 可 靠 性 灵敏 度 的 贡献 来 分 配 ， 并 通过 迭代 模拟 计算 来 得 到 . 本 文 推导 了 基 

于 6 面 截断 重要 抽样 法 的 可 靠 性 灵敏 度 估 计 值 方差 和 变异 系数 的 计算 公式 ， 并 将 该 方法 推广 应 用 

到 并 联系 统 中 . 算 例 结果 表明 :在 估计 值 相 对 误差 小 于 2%、 可 靠 性 灵敏 度 估计 值 变异 系数 相同 

时 ， 基 于 6 面 的 截断 重要 抽样 法 的 可 靠 性 灵敏 度 估计 所 需 的 样本 数 比 传统 重要 抽样 法 、B 球 截断 
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1 引言 


可 靠 性 灵敏 度 常 采用 结构 失效 概率 对 基本 变量 分 布 参数 的 偏 导 数 来 度量 ， 由 于 可 靠 性 灵敏 
度 对 结构 可 靠 性 优化 设计 极为 重要 ， 因 而 得 到 广泛 关注 [1! 均 ， 可 靠 性 灵敏 度 分 析 方 法 可 分 为 近 
似 解析 法 和 数字 模拟 法 .近似 解析 法 一 般 基 于 改进 的 一 次 二 阶 矩 法 对 可 靠 性 作 灵敏 度 分 析 ， 其 
计算 相对 简单 ， 但 该 方法 不 适应 非 线 性 显著 的 极限 状态 方程 情况 . 虽然 数字 模拟 法 比 近似 解析 
法 有 更 为 广泛 的 应 用 范围 ， 但 基于 Monte-Carlo 法 的 可 靠 性 灵敏 度 分 析 对 小 概率 问题 计算 效率 
很 低 ， 这 在 工程 上 难以 接受 . 为 了 加 快 收敛 速度 ， 降 低 计算 量 ， 研 究 人 员 对 Monte-Carlo 法 作 
了 很 多 改进 ， 最 为 常见 的 有 重要 抽样 法 [人 +1、 线 抽样 法 4、B 球 重要 截断 抽样 法 [13, 均 等， 
这 些 方法 在 结构 可 靠 性 灵敏 度 分 析 上 均 有 很 好 的 应 用 ， 其 中 重要 抽样 法 研究 最 为 广泛 . X 
献 [8] 从 提高 投 点 效率 出 发 ， 将 重要 抽样 函数 的 密度 中 心 放 在 极限 状态 方程 的 设计 点 处 ， 从 而 
使 得 在 极限 状态 方程 非 线 性 程度 不 太 大 时 ， 落 入 失效 域 的 样本 点 可 以 达到 50%. 文献 THESE 
了 基于 重要 抽样 法 下 可 靠 性 灵敏 度 估计 量 的 方差 和 变异 系数 ， 给 出 了 可 靠 性 灵敏 度 估 计量 在 给 
定 置信 度 下 的 置信 区 间 . 在 标准 正 态 空间 中 ， 以 原点 为 球 心 ， 可 靠 度 指标 6 为 半径 ， 建 立 8 超 
球面 将 变量 空间 分 割 成 8 超 球 外 区 域 和 6 超 球 内 区 域 ， 落 入 6 超 球 内 的 样本 点 无 需 代 入 功能 
数 计算 其 响应 值 即 可 断定 其 在 安全 域内 ， 这 种 基于 68 超 球 的 截断 抽样 法 在 很 大 程度 上 减少 了 
计算 功能 函数 的 次 数 .文献 [14 综合 重要 抽样 和 6 超 球 截断 抽样 的 优点 ， 提 出 了 一 种 基于 6 球 
的 截断 重要 抽样 法 ， 通 过 引入 以 均值 点 为 球 心 、 可 靠 度 指标 为 半径 的 6 超 球 ， 并 将 重要 抽样 
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区 域 控制 在 8 超 球 外 ， 从 而 减少 了 重要 抽样 函数 在 结构 安全 域 的 抽样 ， 提 高 了 抽样 效率 . X 
献 f1 和 还 将 该 方法 推广 到 多 失效 模式 串联 系统 的 可 靠 性 灵敏 度 分 析 中 ， 推 导 了 可 靠 性 灵敏 度 估 
计 值 的 方差 和 变异 系数 计算 公式 . 与 传统 重要 抽样 法 相 比 ， 该 方法 对 低 维 问题 计算 量 减 少 可 以 
达到 50%. 然而 该 方法 存在 两 个 方面 的 不 足 : 1) 必须 得 到 准确 的 可 靠 度 指标 8 值 ， 如 果 估 算 
的 8 球 半 径 大 于 实际 可 靠 度 指标 ， 会 导致 可 靠 性 灵敏 度 存在 估算 误差 ，2) 研究 对 象 维 数 越 高 ， 
该 方法 效率 越 低 n31， 本 文 针对 上 述 不 足 之 处 ， 提 出 一 种 基于 8 面 分 界 的 重要 抽样 法 ， 来 求解 单 
失效 模式 结构 可 靠 性 灵敏 度 ， 该 方法 在 传统 的 重要 抽样 方法 基础 上 ， 在 设计 点 处 作 关 于 失效 域 
的 切面 ， 即 8 截面 ， 将 变量 空间 分 割 成 重要 抽样 区 域 尺 和 非 重要 抽样 区 域 S， 采 用 基于 C 面 截 
断 的 重要 抽样 函数 抽取 样本 ， 区 域 样本 量 依据 该 区 域 对 可 靠 性 灵敏 度 贡 献 额 来 分 配 ， 而 贡献 额 
在 抽样 过 程 通过 迭代 模拟 得 到 . 算 例 表明 ， 在 相同 的 计算 精度 下 ， 基 于 8 面 的 截断 重要 抽样 法 
所 需 的 样本 数 为 传统 重要 抽样 法 的 50%. 


2 基于 重要 抽样 法 的 可 靠 性 参数 灵敏 度 分 析 


由 于 独立 正 态 变量 具有 一 定 的 普遍 性 ， 因 此 本 文 主要 研究 独立 正 态 变 量 情况 ， 即 假 

设 所 研究 问题 包含 n 维 基本 向 量 x = (zi,zz…… ,Zn) 相 互 独立 且 均 服从 正 态 分 布 ,，x = 
(21,29,:-- En) (i 二 1,2,… ,n), ui fl o2 2: 9WA) x2 3948 S b EZB AA SCIL 

2.1 基于 Monte-Carlo 数字 模拟 法 的 可 靠 性 灵敏 度 分 析 

设 极 限 状 态 变 量 g(x) 是 nn 维基 本 随机 向 量 x 的 函数 ， 则 失效 域 下 和 失效 概率 Py 可 表示 为 


F= {x:g(x) <0), (1) 


P; - P{x:xeF}= | feodx (2) 





其 中 f(x) 是 基本 随机 向 量 x 的 联合 概率 密度 函数 . 由 于 可 靠 性 灵敏 度 可 定义 为 失效 概率 对 基本 
变量 分 布 参数 的 偏 导 数 ， 因 此 (2) 式 对 第 i 个 基本 变量 zi 的 分 布 参数 0; (9; = ni, o.) 的 可 靠 性 灵 
敏 度 如 下 


ðP; = f 2 x); a= f 1709210 jns Ie(x) Of (x) 











x)dx, i= 12, Tha (3) 











00; ðb; ôb; p f(x) 96; 
KeBDOUEEMAEBUSIRB, Ip(x) 为 满足 


T | 0，g(x) >0, 
1, g(x) < 0， 


的 示 性 函数 . 
由 数理 统计 知识 8， 可 推出 (3) 式 基 于 Monte-Carlo 数字 模拟 的 无 偏 估计 值 0 P, /00; 及 其 方 
差 Var(8Py /00;) 分 别 如 下 





La Tr( (xj) ) Of( (x;) 
2» i z1,2,- 4 
^N : f(x) 00, , t , » D, ( ) 








N 
ve( A) e i D r a ho fe) me] 0 
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其 中 xj 为 来 自 于 联合 概率 密度 函数 f(x) 的 第 ;个 子 样 ，N 为 抽样 总 量 . 

2.0 ”基于 重要 抽样 的 可 靠 性 灵敏 度 分 析 

Monte-Carlo 法 收敛 速度 仅 为 O(1/VN)， 当 研究 对 象 为 小 概率 问题 时 ， 需 要 很 大 样本 量 才 
使 (4) 式 得 到 收敛 的 结果 . 为 了 加 快 算法 的 收敛 速度 . 文献 [8] 将 重要 抽样 函数 h(x) 引入 到 可 
靠 性 灵敏 度 分 析 中 ，h(x) 选用 与 基本 随机 向 量 概率 密度 函数 f(x) 相近 的 样式 ， 即 保持 重要 抽 
样 函 数 h(x) 的 方差 和 原 概率 密度 函数 的 方差 一 致 .文献 [10] 研究 表明 ， 重 要 抽样 函数 h(x) 的 
形式 确定 其 抽样 效率 ， 对 于 不 同 的 分 布 参 数 0;， 可 靠 性 灵敏 度 积分 计算 贡献 大 的 样本 点 存 有 差 
异 ， 即 对 于 不 同 分 布 参 数 9; 而 言 ， 可 靠 性 灵敏 度 积分 运算 所 要 求 的 重要 抽样 函数 的 密度 中 心 
是 不 同 的 ， 如 果 对 每 个 分 布 参 数 0; 均 构造 其 相应 的 最 佳 重 要 抽样 函数 ， 这 需要 花费 很 大 的 工 
作 量 来 获得 重要 抽样 法 的 密度 中 心 。 如 果 将 重要 抽样 函数 的 密度 中 心 放 在 max(Ir(x)f(x)) 对 
应 的 样本 ， 即 结构 最 可 能 失效 点 ， 对 每 个 分 布 参数 可 靠 性 灵敏 度 的 积分 计算 而 言 ， 这 样 构 
造 重 要 抽样 函数 的 密度 中 心 不 是 最 优 ， 但 只 要 一 次 寻找 到 max(Ip(x)f(x)) 对 应 的 样本 点 ， 
就 可 建立 重要 抽样 概率 密度 函数 来 计算 所 有 分 布 参数 的 可 靠 性 灵敏 度 . 对 单 失 效 模式 而 
言 ，max(7F(x)F(x)) 所 处 的 样本 ， 对 应 于 极限 状态 方程 的 设计 点 ， 这 可 以 通过 改进 一 次 二 矩 
法 求 得 . 

引入 重要 抽样 密度 函数 h(x)， 则 (3) 式 的 可 靠 性 灵敏 度 可 改写 为 

















ap _ f Ir(x)0f(x) ..[iro990fQ0] .. 
86, = a h(x) 86; h(x)dx — E| h(x) 80; I 2 一 1,2,--- ) 7， (6) 
其 中 重要 抽样 密度 函数 (xc) 满足 
f Noches de 
D 


由 (6) 式 可 以 推导 出 基于 重要 抽样 法 的 可 靠 性 灵敏 度 无 偏 估计 6Py/60; 及 其 方差 Var(BPr/ 
80;) 分 别 如 下 











N 
9P,/00; = 22755) 2 TP (7) 
1 A pel) af)? p1 S Ir) 9f)? 
Var(9Pj /08.) = ni 55. | UNE hoy) 90 |» 6 


J 
其 中 xj 为 来 自 于 联合 概率 密度 函数 f(x) 的 第 ;个 子 样 ，N 为 抽样 总 量 


3 ”基于 6 面 截断 重要 抽样 法 的 可 靠 性 参数 灵敏 度 分析 


3.1 基于 6 面 截断 重要 抽样 法 可 靠 性 灵敏 度 估算 

基于 上 述 研 究 ， 本 文 在 设计 点 Dx = (Dza Dzs,… ,Dz,) 处 作 失 效 面 的 虚拟 切面 ， 
即 B 面 ，8 面 将 变量 空间 刀 分 解 成 重要 区 域 丸和 非 重 要 区 域 5S， 如 图 1 所 示 ， 虚 拟 切 面 的 
含义 是 不 需要 写 出 切面 的 方程 ， 而 直接 由 


R= fx: J (fi 23/0)" - (le: - Dada) - 8 2:0). (9) 
i i-i 


i=1 
S= pe: 35 (s nae)" - 3, (e= Diod =p «tj (10) 


i= 
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所 示 的 几何 关系 来 确定 玉 和 3 区 域 ， 这 对 于 隐 式 函数 来 说 可 以 减少 较 大 的 计算 工作 量 ， 因 为 作 
切面 需要 求 隐 函数 的 偏 导数 . 


g(X)20 


[i planc 








图 1: 重要 抽样 区 域 忆 与 非 重要 抽样 区 域 3 示 意图 


将 变量 空间 由 虚拟 切面 6 面 划分 为 尺 和 53 区域 后 ， 可 靠 性 灵敏 度 9Pr/6b 可 由 整个 变量 空 
闻 呈 的 失效 概率 分 解 成 及 和 3 区 域 的 失效 概率 9Ppyj/66 M 0 Ps ;/00; 之 和 . 


0P,/00; = Prf /00; + OPs,/00; 











Ir(x) Of (x) Ip(x) 9f (x) 
n ha(x) 86; hg(x)dx E hs(x) 80; h (x)dx, (11) 
其 中 Ag (x) M hs (x) 为 截断 重要 抽样 函数 ， 可 取 为 以 下 形式 

0, xéR, 
hr(x) = | á (12) 

h(x)/f,h(x)dx, x€ R, 
x)= " Xs (13) 

h(x)/f;h(x)dx, x€ S, 


故 可 靠 性 灵敏 度 基于 9 面 截断 重要 抽样 的 无 偏 估计 为 
9P,/00; = OPrs/00; + OPss/00; 








N Ns 
R 1 Ir(xR ) of Om) g 1 Ip (xs;) Of (xs,) , 
一 —1,2,---,n, (14 
QD S) S + 2- Ness) 96; i moe 


其 中 xm xs; 分 别 来 自 8 面 截断 重要 抽样 函数 h(x), hs(x) 的 第 j 个 样本 , Ne, Ns 2I n9), 
hs(x) 的 样本 数量 .在 独立 正 态 变量 情况 下 ，n 维 变量 的 联合 概率 密度 函数 f(x) 可 改写 成 每 个 
变量 边缘 概率 密度 函数 f) (i = 1,2, ,n) 的 乘积 形式 ， 即 


= ]I52. 


180 I 8 X 学 学 R 第 28 卷 











H.0; H5 fila) AR. 当 0;= ji 时 











Oft) IF _ Bhi o gala. 
00; Opi c2 py chou (15) 
= 9; — coi 时 ， 
afe) Of)  lp(m-nmy oo 
00; z 0c; e gi [( gi ) 1] f(x), i= 1,2, QD, (16) 


将 (15), (16) 式 代入 (14) 式 ， 得 到 可 靠 性 灵敏 度 估计 值 9Pr/8w M OP, / 0o; 的 具体 表达 式 如 下 


Ne 

^ B f (xR;) Ir(xn;) zi — Hi 

dE. X Na hn(xm) o? 
ge 








y X f(xs;) Ir(xs;) zi — ui 


, i—12,-.,n, 17 
Ns hs(xs) o? qn) 


Sj;=1 








ES ea oq 
Rj-i 


Ng hg(xnj)ci gi 





N. 
和 f(xs;) Ir(xs;) 1 p yzi — piy? 
ex Ns cr Oi ) il. (18) 
3.0 基于 0 面 截 断 重要 抽样 可 靠 性 灵敏 度 估计 值 的 方差 分 析 
采用 数字 模拟 的 方法 对 真 值 进行 估计 均 是 近似 的 ， 为 了 对 基于 6 面 截断 重要 抽样 法 的 收敛 
性 有 所 了 解 ， 有 必要 对 估计 值 作 方差 分 析 . 对 (11) 式 所 示 的 估计 值 求 数学 期 望 和 方差 可 得 

















som 
Var [3 P; /00;| = Var ES ME | 十 Var P = 2r) | ; (20) 


在 数值 模拟 的 过 程 中 考虑 用 样本 平均 值 和 方差 代替 总 体 的 数学 期 望 和 方差 ， 可 近似 得 到 可 靠 性 
灵敏 度 估 计 值 9Py/69; 的 方差 和 变异 系数 分 别 如 下 


Var(OP, /80;) = Var(OPn,/00) + Var(0Ps; /00:) 


Ng 


1 «A (Ir(xn) Of(xn;) i? Le pd) 2 
xen P ever 86; xl » R(xn,) ð; 2 md 














1 SS p Ir(xs;) 0f(xs;) i2 «3 Ir(xs, 1/00 
+ Ns(Ns — 1) P» ð; A? 2. , hs(xs;) } |: (21) 


Cov(OP,/80;)— 4/ Var(3Êş /00:) / [oP;/00j, i—1,2,-,n. (22) 
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3.3 基 0 面 截断 重要 抽样 的 实施 

以 下 将 讨论 8 面 截 断 重要 抽样 法 实现 过 程 中 的 两 个 关键 问题 ， 分 别 是 8 面 截断 重要 函数 抽 
样 问题 和 区 域 、 合 适 的 样本 量 问题 . 

1) B 面 截断 重要 函数 抽样 问题 

采用 h(x) 抽取 样本 点 x， 如 果 


n 


> (Gi = pi)/0i) - M5 (i - D4/o)? - £ 20 
i-l i=l 

WE. WUFRARGXSENGESEDONDR, TU. $ ri = 2D。 一 zi (i = 1,2,… ,n)， 新 得 到 的 x 即 落 
TEE EDOBR R, 即 为 hrn(x) 的 样本 ， 具 体 见 图 2 中 左 图 所 示 所 示 . 同 理 ， 对 于 非 重要 区 域 5， 
采用 js(x) 抽样 抽取 样本 点 x， 如 果 


T, 


(Gi mi)/oi) — Y (Gi Du)/oi) -8 «0 

i-l1l i 
满足 ， 则 样本 点 落 入 非 重 要 区 域 S， 即 为 hs(Co FERES. BW, Ae = 2D, -rli = 
1 2,…… ,n)， 新 得 到 的 x 即 落 在 非 重 要 区 域 S， 即 为 hs(x) 的 样本 ， 详 见 图 2 中 的 右 图 所 示 . 


i 
ra 


4 4 
AR © 























图 2: DB 截面 截断 重要 抽样 函数 抽样 示意 图 


2) ERR, 8 合适 的 样本 量 问题 
R, S 区 域 合适 的 样本 量 ， 可 采用 迭代 抽样 来 获取 ， 具 体 过 程 可 表述 如 下 ， 
当 抽样 总 数 N 确 定时， 文献 115 研究 表明 ， 按 VVar(Ppy) 与 / Var(Psy) 的 比率 来 分 配子 
空间 RR 与 5 的 抽样 量 Na 和 Ns， 可 使 估计 值 方差 最 小 ， 因 此 可 按 








Nr=N Var(Prs) / (VVar(Pas) 十 V Var(Ps;) J (23) 





Ns = NJ Va(Ps p) / (VVar(Pas) + V Var(Ps;) 上 (24) 


来 确定 Ng 和 Ns. 
HF Var(Êrs), Var(Êsp 均 为 未 知 ， 因 此 本 文采 用 和 迭代 抽样 方式 ， 通 过 已 有 抽样 信息 ， 来 
确定 下 一 次 R, 5S 区 域 合适 抽样 量 ， 具 体 过 程 如 下 : 
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设 前 1 次 在 R 和 5 区 域 的 样本 量 分 别 为 NO 和 NN 人 ) (I = 0,1, ,有 JU Var( P0), Var( PO) 
可 以 表示 为 














DONE f(x )TF (Xn) .2 
Var(P (0) — : 3 | Y { ní RÍ } 
NO NO 4 k=1 p(k)_ hn(xgmo) 
(Z Nà (ÈE Ng )-3) R} =1 à 
k=l k=1 
1 l NE? FRH Mr (xpo) W] (25) 
(X NO) E aea hrx po) 
k=1 
: DONO f(xso)Hr(xgi) 2 
Var(P£7) = E CFL od) jl A hs(Xg%) j 
(22 Ns WCE N$)-1) isa ' 


l NS f(x o» )TF (x (x) ) 


5j sí 2 
| En H >, | js(xsto) jj | (26) 





第 (二 1) 次 抽样 在 及 和 3 区 域 的 样本 量 NOD, NUT? (1 > 0) 可 采用 下 式 确定 
Nf» NOD var) / (| Var( ÊH) + y varc) )， (27) 


NOD = NO? Var (BY) / (VVar( Po) 十 V var) ), (28) 


其 中 NC 为 第 (! + 1) AARETE [i] fl RH FE i. 

3.4 基于 0 截面 的 截断 重要 抽样 法 的 估算 步骤 

根据 上 述 研究 ， 可 以 给 出 基于 6 截面 重要 抽样 法 求解 可 靠 性 灵敏 度 的 具体 实施 步骤 如 下 : 

1) 通过 一 次 可 靠 度 方法 求解 得 到 单个 极限 状态 方程 的 设计 点 De 和 可 靠 度 指标 6; 

2) 在 设计 点 Dx 处 作 失 效 面 的 虚拟 切面 ， 即 8 面 ， 将 变量 空间 分 割 成 重要 积分 区 域 RUE 
重要 抽样 区 域 $， 并 建立 截断 重要 抽样 函数 hr(x) M hs (x); 

3) 第 一 次 抽样 

假定 Var(Pg)) = Var( BÉ) (上 标 (0) 表示 第 0 次 迭代 运算 )， 通 过 (27), (28) 式 计算 得 到 第 一 
次 采用 截断 重要 抽样 函数 (xc) 和 hs(x) 抽取 的 样本 量 NC. NP, 

4) 第 (1 十 1) 次 抽样 

通过 前 1 次 抽样 信息 ， 代 入 (25), (26) 式 中 计算 Var( ÊG), Var(P82), W (27), (28) 式 计算 
得 到 第 (! + 1) 次 采用 截断 重要 抽样 函数 ha(x) 和 hs(x) 抽取 的 样本 量 N(D. NO, 

5) HI (17), (18) 式 计算 可 靠 性 灵敏 度 的 估计 值 ， 通 过 (21), (22) 式 计算 其 方差 和 变异 系 
数 ， 如 果 变 异 系数 小 于 给 定 值 ， 即 算法 收 伍 ， 否 则 令 1 = 1 十 1， 返 回 第 4) 步 至 算法 收敛 . 








4 算 例 分 析 


算 例 1 高 维 非 线性 极限 状态 函数 为 g(x) = 7z? 十 2 十 2 十 X32 十 2 一 ze 十 3， 各 随机 变量 相 
互 独立 ， 且 服从 标准 正 态 分 布 . 采用 改进 的 一 次 二 阶 矩 法 求 得 设计 点 为 : [0,0,0,0,0,3]， 采 用 
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传统 重要 抽样 法 、B 球 截断 重要 抽样 法 和 6 面 截断 重要 抽样 法 三 种 方法 对 重要 变量 ze 作 可 靠 性 
灵敏 度 分 析 ， 其 可 靠 性 灵敏 度 估计 和 变异 系数 计算 结果 见 表 1. 


Ar 算 例 1 的 计算 结果 














方法 传统 重要 抽样 法 [9 球 截断 重要 抽样 法 有 面 截断 重要 抽样 法 

抽样 数目 5 x 106 4.05x106 5x106  2.65x109 5x109 

估计 值 0.000026 0.000026 0.000026 0.000026 0.000026 

OP;/Oug 变异 系数 0.00642 0.00639 0.00575 0.00635 0.00454 
HITI É (96) / 0.07736 0.2830 1.0838 0.09498 

估计 值 0.000094 0.000093 0.000093 0.000083 0.000094 

OPy/Boe 变异 系数 0.00586 0.00583 0.00525 0.00578 0.00414 
相对 误差 (%) / 0.07282 0.2594 0.9943 0.1063 








以 传统 重要 抽样 法 估算 值 为 精确 值 ， 算 例 1 采 用 本 文 方法 估算 可 靠 性 灵敏 度 ， 其 相对 误 
差 不 超过 2%， 在 估计 值 相同 变 异 系 数 下 ， 本 文 方法 所 需 样本 量 是 传统 重要 抽样 法 的 52%， 
是 8 球 截断 重要 抽样 法 的 64.197%. 

算 例 2 ”如 图 3 所 示 的 九 盒 段 机 辟 模 拟 结构 由 64 个 杆 元 件 和 42 个 板 元 件 构成 ， 材 料 为 铝 合 
金 ， 原 始 数据 取 自 文献 [16]. 各 变量 服从 正 态 分 布 ， 外 载荷 己 的 均值 和 变异 系数 分 别 为 WP = 
150kg, Cov(P) = 0.25， 第 i 个 单元 强度 Ri 的 均值 和 变异 系数 分 别 为 Ri = 83.5 kg, Cov(Ri) = 
0.12. 由 失效 模式 的 枚 举 方法 可 求 得 结构 体系 的 一 个 主要 失效 模式 的 极限 状态 函数 为 g 
4.0R6s — 3.9998Rz7 十 4.0R7s 一 .可靠 性 灵敏 度 估 计 值 及 其 变异 系数 如 表 2 所 示 ， 

以 传统 重要 抽样 法 估算 值 为 精确 值 ， 算 例 2 采 用 本 文 方法 估算 可 靠 性 灵敏 度 ， 其 相对 
误差 不 超过 1%， 在 估计 值 相同 变异 系数 下 ， 算 例 2 本文 方 法 所 需 样本 量 是 传统 重要 抽样 法 
的 63%， 是 8 球 截 断 重要 抽样 法 的 8496. 






180cm 


"t (0 2pM'^ — 2R/. 
20cmT?9 人 


[MM ERE 





图 3: 九 盒 段 结构 示意 图 


算 例 3 ”并 联系 统 由 两 极限 状态 方程 组 成 gl(x) = z? — 52i 一 82 - 16, g(x) = 一 16z1 + 
zZ2 十 32， 其 中 包含 的 随机 变量 均 服 从 相互 独立 的 标准 正 态 分 布 ， 先 通过 优化 方法 得 到 并 联系 统 
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最 可 能 失效 点 ， 以 该 点 为 抽样 中 心 ， 按 照 (12), (013) 式 建立 截面 截断 重要 抽样 函数 ， 其 可 靠 性 
灵敏 度 估算 见 下 表 3 Pp. 

算 例 3 显示 本 文 所 提 方 法 可 以 推广 应 用 到 并 联系 统 可 靠 性 灵敏 度 分 析 中 ， 与 传统 重要 抽样 
法 相 比 ， 本 文 方法 估算 失效 概率 的 相对 误差 小 于 1%， 在 估计 值 相 同 变异 系数 下 ， 本 文 方法 所 
需 计算 量 为 传统 重要 抽样 法 的 46%， 是 8 球 截断 重要 抽样 法 的 77.966%. 


表 2: 算 例 2 的 计算 结果 
































方法 传统 重要 抽样 法 BRANNER k J 面 截断 重要 抽样 法 

抽样 数目 5 x 109 3.755 x100 5x10 3135x106 5x109 

估计 值 -0.00134 -0.00133 -0.00134 -0.0013 -0.00134 

OP;/Ouss 变异 系数 0.00089 0.00089 0.00077 0.00094 0.00075 
相对 误差 (%) / 0.01572 0.07413 0.02396 0.03220 

估计 值 0.00158 0.00158 0.00158 0.00158 0.00158 

OP;/0css 变异 系数 0.00149 0.00149 0.00129 0.00169 0.00134 
相对 误差 (%) / 0.02781 0.1157 0.04297 0.03033 

估计 值 0.00133 0.00133 0.00133 0.00133 0.00133 
OPj/üwr 变异 系数 0.00087 0.00087 0.00075 0.00090 0.00071 
相对 误差 (96) / 0.01501 0.02551 0.02176 0.01200 

估计 值 0.00158 0.00158 0.00158 0.00158 0.00158 

8 户 /6crr 变异 系数 0.00144 0.00144 0.00124 0.00159 0.00126 
相对 误差 (%) / 0.00807 0.06721 0.09257 0.06721 

估计 值 -0.00134 -0.00134 -0.00134  -0.00133 -0.00133 

OP;/Onrs 变异 系数 0.00084 0.00083 0.00072 0.00081 0.00065 
相对 误差 (96) / 0.00974 0.00749 0.09811 0.06816 

估计 值 0.00158 0.00158 0.00158 0.00158 0.00158 

OP;/0cs 变异 系数 0.00136 0.00135 0.00117 — 0.00144 0.00144 
相对 误差 (96) / 0.01011 0.01517 . 0.2415 0.1523 

估计 值 0.00334 0.00334 0.003334 — 0.00334 0.00334 

OP;/Oup 变异 系数 0.00095 0.00094 0.00082 0.0011 0.00083 
相对 误差 (号 ) 0.01499 0.04227 . 0.06415 0.01001 

估计 值 0.00037 0.00037 0.00037 0.00037 0.00037 

OP;/0cop 变异 系数 0.00169 0.00168 0.00146 — 0.00120 0.00159 
相对 误差 (96) / 0.01920 0.08464 0.07436 0.1444 











本 文 在 传统 重要 抽样 法 基础 上 提出 了 基于 6 面 的 截断 重要 抽样 法 来 求解 结构 可 靠 性 灵敏 
度 .该 方法 在 设计 点 处 作 失 效 面 的 虚拟 切面 ， 即 68 截面， 将 变量 空间 分 成 重要 抽样 域 尺 和 非 重 
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要 抽样 域 9S， 并 在 丸和 3 区 域 分 别 建立 相应 的 截断 重要 抽样 函数 ， 根 据 丸 和 3 区 域 对 可 靠 性 灵 
敏 度 的 贡献 额 不 同 来 分 配 截面 截断 重要 抽样 函数 的 样本 量 ， 从 而 提高 了 抽样 效率 . 算 例 分 析 表 
明 ， 与 重要 抽样 法 估算 结果 相 比 ， 本 文 方法 的 相对 误差 小 于 2%;， 在 估计 值 相同 变异 系数 下 ， 
该 方法 毛 需 样本 量 比 传统 重要 抽样 法 少 近 50%; 与 B 球 截断 重要 抽样 法 相 比 ， 本 文 方法 在 高 维 
情况 下 有 明显 优势 ， 并 且 并 不 需要 准确 求 得 可 靠 度 指 标 ， 因 而 本 文 方法 应 用 范围 更 广 ， 本文 方 
法 还 可 以 推广 应 用 到 并 联系 统 可 靠 性 灵敏 度 分 析 中 . 


A3: 算 例 3 的 计算 结果 

















方 法 传统 重要 抽样 法 有 B 球 截断 重要 抽样 法 DB 面 截断 重要 抽样 法 

抽样 数目 5 x 108 2.95 x 10 5x10 23x106 5x10? 

估计 值 0.00398 0.00398 0.00398 0.00398 0.00398 

OP;/0n| 变异 系数 0.00134 0.00133 0.00102 0.00132 0.00089 
相对 误差 (%) / 0.01308 0.03069 0.1185 0.06037 

估计 和 值 0.00857 0.00857 0.00858 0.00858 0.00858 

OP;/0c 变异 系数 0.00122 0.00122 0.000933 0.00119 0.00081 
73t E (90) / 0.01423 0.03826 ^ 0.1015 0.06334 

45 THÉ 0.00257 0.00257 0.00257 0.00257 0.00257 

OP,/Ouo 变异 系数 0.00131 0.00131 0.00100 — 0.00130 0.00088 
J8x] iX E (96) / 0.01208 0.02455 0.1122 0.05339 

估计 值 0.00277 0.00277 0.00277 0.00277 0.00277 

Ə; /ðo 变异 系数 0.00120 0.00119 0.00091 — 0.00116 0.00079 
相对 误差 (%) / 0.01119 0.01553 0.08016 0.04008 
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Improved Reliability Sensitivity Estimation and its Variance Analysis 
by a Novel 8 Hyper-plane Based Importance Sampling Method 


ZHANG Feng?, LV Zhen-zhou?, CUI Li-jie? 


(1- College of Mechanics and Civil & Architecture, Northwestern 
Polytechnical University, Xi'an 710129; 2- College of Aeronautics, Northwestern 
Polytechnical University, Xi'an 710072) 


Abstract: A novel 8 hyper-plane based importance sampling method is presented to estimate relia- 
bility sensitivity of a structure. By introducing a virtual hyper-plane tangent to the failure surface, the 
variable space is separated into an importance region R and a unimportance region S, on which the 
truncated importance sampling functions hg(x) and hs(x) are established, respectively. The sampling 
numbers generated from hn(x) and hs(x) are dependent on the contribution of the reliability sensitivity, 
which is determined by the iterative simulations. The formulae of the reliability sensitivity estimation, 
the variance and the coefficient of variation are derived for the presented 8 hyper-plane importance 
sampling method. The presented method is suitable for the reliability sensitivity estimation of both 
the single failure mode and the multiple failure mode in parallel. Examples show that the proposed 
method is more efficient than the traditional importance sampling method and the 8 hyper-sphere 
importance sampling method, in the case that the variation coefficients of three estimations keep the 
same quantity and the relative errors of the reliability sensitivity estimations are less than 2%. 
Keywords: design point; 8 hyper-plane; 9 hyper-plane based importance sampling method; reliability 
sensitivity; variation coefficient 
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